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(54) Title: METHOD FOR LOW-FUEL. COMPUTER-ASSISTED CONTROL OF ANY NUMBER OF THRUSTERS ARRANGED ON 
A SPACECRAFT 

(54) Bezeichnung: VERFAHREN ZUR TREIBSTOFFMINIMALEN, RECHNERGESTOTZTEN ANSTEUERUNG FOR BELIEBIG AN 
EINEM RAUMFAHRZEUG ANGEORDNETE DUSEN 

(57) Abstract 

The invention relates to a method for low-fuel, computer-assisted control of any number of thrusts mounted in any given arrangement 
on a spacecraft, which solves the problem of linear optimization by reducing computational effort in a very short time. The method is based 
on the dual simplex algorithm and suitable provision of a dual admission start table. The dual admission start table is set up on the basis of 
a previously calculated optimum table and a transformation of the actual force/power vector pertaining to the chosen optimum table. Given 
processing steps of the dual simplex algorithm are skipped totally or are substituted by processing steps which require substantially less 
computational effort. Finally, a reduced table as opposed to the simplex table is used. Said reduced table consists solely of an identification 
number for selecting a thruster from the entire number of thrusters and the vector of the control values of the thnisters which, arc contained 
in said selection. Reduction in computational effort is achieved by carrying out given processing steps in advance and by providing the 
corresponding results in the form of data fields and index or identification number matrixes in the working memory of the computer. Said 
results torn the initialization phase of the method need only to be addressed and retrieved in order to calculate thruster control. 

(57) Zusanmienfassung 

Es soli ein Verfahren zur insgesamt treibstoffminimalen, rechnergestiltzten Ansteuerung einer beliebig festgelegten Anzahl und beliebig 
festgelegten Anordnung von Dasen an einem Raumfahizeug angegeben werden. das dieses lineare Optimierungsproblem durch Reduzierung 
des Rechenaufwands in maglichst kurzer Zeit 16st. Das Verfahren basiert auf dem dualen Simplexalgorithmus und der geeigneten 
Bereitstellung eines dual zulflssigen Starttableaus. Das dual zulSssige Starttableau wird aus einem vorab bercchneten Optimaltableau 
und der zu diesem gewShlten Starttableau dazugeh5rigen Transformation des aktuellen Ki^fte-Momenten-Vektors aufgebaut Bestimmte 
Verarbeitungssschritte des dualen Simplexalgorithmus entfallen vollst^dig bzw. werden durch Verarbeitungsschritte substituiert, die 
wesentlich weniger Rechenaufwand benOtigen. SchlieBlich wird nicht mit dem Simplextableau, sondem mit einem reduzierten Tableau 
gearbeitet. Dieses reduzierte Tableau besteht lediglich aus einer Identifikationsnummer fUr eine Dtlsenauswahl aus der Gesamtzahl der Dilsen 
und dem Vektor der Ansteuerwerte flir diejenigen Dtlsen, die in dieser DOsenauswahl enthalten sind. Die Reduzierung des Rechenaufwands 
wird dadurch erreicht, daS bestimmte Verarbeitungsschritte vorab durchgefQhrt werden und die entsprechenden Ergebnisse dann in Form 
von Datenfeldern und Indizes- bzw. Identifikationsnummemmatrizen im Arbeitsspeicher des Rechners bereitliegen. Diese Ergebnisse des 
Initialisierungstesis des Verfahrens brauchen dann zur Berechnung der DQsenansteuerung nur noch adressiert und abgerufen zu werden. 
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Verfah ren zur treibs tof f minimalen , rechnerqesttitz ten Ansteue- 
rung far bel iebiq an einem Raumfahrzeuq anqeordnete DUsen 



15 Die Erfindung bezieht sich auf ein Verfahren zur insgesamt 

treibstof fminimalen, rechnergestiitzten Ansteuerung einer 
beliebig festgelegten Anzahl und beliebig festgelegten 
Anordnung von DUsen an einem Raumf ahrzeug. Die Diisen konnen 
dabei diskret Oder auch kontinuierlich ansteuerbar sein. 

20 

Fiir die Ausubung von auOeren Kraften und Momenten an einem 
Raumfahrzeug ist der Einsatz von Diisen in vielen Situationen 
notwendig. Da der Einsatz der Dusen Treibstoff verbraucht und 
dieser unmittelbar Einfluli auf das Startgewicht und die 

25 Lebensdauer von Raumf ahrzeugen und insbesondere von Satelliten 

hat, koromt der DUsenans teuerung mit in der Summe minimalem 
Treibstof fverbrauch eine wichtige Bedeutung zu. DUsenanordnun- 
gen, die prinzipiell einen niedrigen Treibstof fverbrauch 
ermoglichen, eine hohere Oder gleiche Redundanz bei kleinerer 

30 DUsenanzahl aufweisen und/oder durch die Nutzlast am Raumfahr- 

zeug beeinfluBt sind, sind oft unregelmMBige und schiefwinklige 
Dtisenanordnungen. Aus diesen Grunden ist die Forderung nach der 
treibstof fminimalen Diisenansteuerung, insbesondere auch filr 
beliebige (schiefwinklige) feste Diisenanordnungen an Raumf ahr- 

35 zeugen von besonderer Relevanz. 

Es ist bekannt, daB sich das Problem der treibstof fminimalen 
DUsenansteuerung auf ein lineares Optimierungsproblem zurttck- 
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fUhren IfiBt, fiir das der Slmplexalgori thmus ein geelgnetes und 
effizientes Ldsungsverfahren darstellt (siehe zum Beispiel 
P.J. Wiktor **Minlmum Control Authority Plot: A Tool for 
Designing Thruster Systems" In: Journal of Guidance, Control, 
and Dynamics, Vol* 17, No. 5, Sept. -Oct. 1994, pp. 4-5). Der 
Einsatz des Simplexalgorithmus scheitert jedoch haufig an der 
ben&tigten Rechenzeit. So kdnnen oft die Echtzeitanf orderungen 
zum Beispiel fiir den Einsatz des Simplexalgorithmus innerhalb 
des Regeltaktes des Lage- und Pes it ions regelkreises eines 
Raumf ahrzeuges nicht erfiillt werden, well die besonders in 
hohen Orbits einsetzbaren Rechner den Simplexalgorithmus nicht 
ausreichend schnell abarbeiten konnen. 

Eine Vielzahl der sich heute im Einsatz befindenden Raumfahr- 
zeuge (Satelliten) besitzen eine Diisenanordnung, bei der jede 
DUse zum korperfesten koordinatensystem des Raumf ahrzeugs zu 
einer Achse parallel und damit zu den zwei anderen Achsen 
orthogonal angeordnet ist. Diese orthogonale und 
achsenparallele Anordnung hat zur Folge, daB jede DUse nur 
Krfifte in einem Richtungssinn einer Achsenrichtungr und Momente 
nur in einem Drehsinn urn die beiden anderen Achsen erzeugen , 
kann. Der Vorteil dieser Anordnung ist, daB die Steuerung eines 
jeden Preiheitsgrades des Raumf ahrzeuges als starrem Korper 
durch den Einsatz von jeweils nur zwei Diisen moglich ist, ohne 
andere Freiheitsgrade zu beeinf lussen. Die Ansteuerung der 
Diisen ist fiir diese Diisenanordnung sehr einfach. Der Nachteil 
dieser Diisenanordnung besteht darin, daB mindestens zwolf Diisen 
notwendig sind, um alle Freiheitsgrade gleichzeitig und 
beliebig steuern zu konnen. MuB der Ausfall von zwei beliebigen 
Diisen tolerierbar sein, sind 36 Diisen erf orderlich. 

Demgegeniiber liegt der Vorteil einer geeigneten schiefwink- 
ligen, also nicht orthogonalen und achsenparallelen Diisen- 
anordnung Im wesentlichen in folgenden drei Punkten: 

1. Bei frei wShlbarer fester Diisenanordnung sind weniger als 

12 Diisen erf orderlich. 

2. Fiir die Tolerierung des Ausfalls von zwei beliebigen Diisen 
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sind erheblich weniger als 36 Dusen erf orderlich. 
3. Bei trelbstoffmlnimaler Ansteuerung der DUsen Ist eln 
geringerer Tre lbs toff verbrauch moglich. 

Der Nachteil dieser Dusenanordnung besteht darin, dal^ die 
treibstof f minimale Dlisenans teuerung die Losung eines linearen 
Optimierungsproblems erfordert und der Rechenauf wand ftir die 
Losung dieses Problems mit dem Simplexalgorithmus oft zu hoch 
ist. Aus diesem Grunde sind andere Verfahren entwickelt worden, 
die weniger Rechenzeit benotigen und durch die eine mehr Oder 
weniger treibstof fminimierte , aber keine insgesamt treibstoff- 
minimale Dusenansteuerung erreicht werden kann* 

Beim sogenannten **Linear Programming'* wird das Diisenansteue- 
rungsproblem als lineares Optimierungsproblem aufgefaBt und mit 
einschlagigen Verfahren der linearen Optimierung, wie sum 
Beispiel den Simplexalgorithmus , gelost. Der schon erwMhnte 
Nachteil, der den Einsatz dieses Verfahrens in vielen FMllen 
unmoglich macht, ist der zu hohe Rechenauf wand. 

Beim ••Table Look-up" -Verfahren werden alle treibstof fminimaLen 
Diisenanordnungen fvir Einheitskraf te und -momente bezUglich 
einer Achse des korperfesten Koordinatensys terns des Raumfahr- 
zeugs berechnet (Losen des linearen Optimierungsproblems) und 
in einer Tabelle abgelegt. Sind im Einsatz Krafte und Momente 
bezuglich mehrerer Achsen gleichzeitig auszuxiben, werden die 
entsprechenden Dusenansteuerungen bezUglich Kraft Oder Moment 
und fiir jede geforderte Achse und dem geforderten Richtungs- 
Oder Drehsinn aus der Tabelle ausgelesen und mit dem entspre- 
chenden Betragswert der Kraft- bzw. Momentenkomponente in 
dieser Richtung multipliziert • Den Ansteuerwert ftir jede DUse 
erhalt man durch die Addition der einzelnen Werte ftir die 
jeweilige Dtise. Dieses Verfahren ist sehr einfach und sehr 
schnell, ftihrt im allgemeinen aber zu nicht treibstof fminimalen 
Dttsenansteuerungen. 

Bei veffeinerten und erweiterten ••Table Look-rup'* -Verfahren 
handelt es sich um Verfahren, die eine Look-up-Tabelle anlegen, 

3 
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die EintrSge fUr zusatzlich eingefUhrte Achsen enthalten. Durch 
diese Verfahren kttnnen treibstoffsparende Diisenans teuerungen 
ermittelt werden. Zugleich nimmt aber die Obersichtlichkeit des 
Verfahrens ab und wegen der zunehmenden Gr&Be der Tabellen der 
5 Rechenaufwand erheblich zu. 

Andere Verfahren beruhen au£ der Auswertung und Verbesserung 
des Ergebnisses, das mit dem "Table Look-up" -Verfahren 
ermittelt wird- Das Ziel dieser Auswertung besteht darin, 

10 besonders treibs tof fverschwendende Diisenans teuerungen zu 

verhindern bzw. geeignet zu korrigieren, Hier soil exemplarisch 
ein Verfahren, das auf einem Fuzzy-Ansatz beruht, genannt 
werden (vgl. T.Suzuki, K.Yasuda, S.Yoshikawa, K.Yamada, 
N.Yoshida "An Application of Fuzzy Algorithm to Thruster 

15 Control Systems of Spacecraft" In: ISTS 94-C-16, 19th Inter- 

national Symposium on Space Technology and Science, Yokohama, 
Japan, May 15-24, 1994, pp. 1 - 8). 

Eine Treibstof f minimierung der Dusenansteuerung insgesamt hat 
20 zum Ziel eine solche Ansteuerung der Diisen zu ermitteln, die 

* 

1* die geforderten Krafte und Momente am Raumfahrzeug ausUbt, 
2. fiir jede Diise einen nichtnegativen Ansteuerwert liefert, 
und 

25 3. die Summe aller Ans teuerwerte der Diisen minimiert. 

Dies stellt - wie erwahnt - ein lineares Optimierungsproblem 
dar. Da das Losen des Problems in. der Kegel zu zeitauf wendig 
ist, kann die Losung des quadratischen Optimierungsproblems als 
30 Ersatzproblem durchgef iihrt werden. 

Die Losung des Ersatzproblems fiihrt zur Ermittlung derjenigen 
Diisenans teuerung, die 

35 1« die geforderten KrMfte und Momente am Raumfahrzeug ausiibt, 

2. fiir jede Diise einen nichtnegativen Ansteuerwert liefert, 
und 



BNSOOCID:<WO fiatSOSAAl I > 



wo 98/49058 



PCT/EP98/01595 



3. die Summe aller quadratlschen Ansteuerwerte der DUsen 
mlnimiert • 

Dieses quadratische Optimierungsproblem kann liber die Bereit- 
5 stellung der Moore-Penrose-invertierten Sys temmatrix der 

Diisenanordnung am Raumf ahrzeug und eines Of f setvektors , der mit 
Hilfe der singulMren Wertezerlegung berechnet wird (vgl. z.B. 
B.Noble^ I^W. Daniel "Applied linear algebra" Prentice-Hall 
International, INC. (1988), pp. 338 - 350), gelost werden. Die 

10 Losung des quadratlschen Optimierungsproblems kommt dann als 

Produkt der Moor*Penrose-invertierten Systemmatrix mit dem 
kommandierten KrM£te-Momenten-Vektor und der Addition des 
0££setvektors zustande. Die Addition des Offsetvektors bewirkt, 
daB.im LSsungsvektor nur nicht negative Ansteuerwerte 

15 auftreten. Die Berechnung des Offsetvektors wird jedoch dann 

numerisch immer aufwendiger, wenn die Zahl der anzusteuernden 
Diisen um zwei oder mehr iiber der Dimension des kommandierten 
Krfif te-Momenten-Vektors liegt. Im allgemeinen ist die Losung 
nicht treibstof f-minimal , fallt aber mit der Losung des 

20 linearen Optimierungs-problems zusammen, wenn die Zahl der 

anzusteuernden Diisen nur um eins groBer ist als die Dimension 
des kommandierten Kraf te-Momenten-Vektors {vergleiche. z.B. 

DE 195 19 732 A1 ) . Weiterhin ist nachteilig, daB langere 
Einschaltzeiten iiberproportional "bestraft" werden und daher 

25 . seltener in den Ans teuerungen vorkommen. Fiir den Betrieb der 

Diisen und ihrem Wlrkungsgrad ist es jedoch giinstiger, sie mit 
IMngeren Einschaltzeiten zu betreiben. 

Der Nachte.il bestehender Verfahren besteht mithin darin, daB 
30 sie entweder nicht das lineare Optimierungsproblem losen und 

damit im allgemeinen nicht die treibstof f minimale Ansteuerung 
der DUsen bereitstellen oder aber, wenn bei den Ansteuerver- 
fahren das lineare Optimierungsproblem gel5st wird und der 
Simplexalgorithmus ausgefUhrt wird, dieses ublicherweise zu 
35 rechenaufwendig ist. 

Der Erfindung liegt die Aufgabe zugrunde, ein Verfahren der 
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eingangs angegebenen Art bereitzustellen. mit dem das lineare 
Optimierungsproblem filr die treibstof f roinimale Ansteuerung 
beliebig fester DUsenanordnungen aur gleichzeltlgen Reali- 
sierung von KrSften und Momenten an elnem Raumfahrzeug in 
signifikant kUrzerer 2eit gelSst werden kann. als es bisher 
mdglich ist. 

Diese Aufg^be wird gemSO der Erfindung durch die im Anspruch 1 
gekennzeichneten Merkmale gel6st. 

Das Verfahren nach der Erfindung basiert also vorteilhaft auf 
einem, dem Simplexalgorithmus angelehnten Verfahren, das jedoch 
aufgrund der konsequenten Ausnutzung der Besonderheiten bei der 
DUsenansteuerung bei Raurof ahrzeugen und der Aufgabe der Verwen- 
dung des Ublichen Simplextableaus zu einer RechenzeitverkUrzung 
von ca. einer GrdDenordnung fUhrt. 

Vorteilhafte Ausges taltungen und Weiterf Uhrungen des Verfahrens 
nach der Erfindung sind in den dem Anspruch 1 nachf olgenden Pa- 
tentansprUchen gekennzeichnet . 

Das Verfahren nach der Erfindung basiert auf folgenden Grund- 
lagen: 

a) Das bei der DUsenansteuerung zu losende lineare Optimie- 
rungsproblem kann wie folgt auf geschrieben werden: 
Gesucht ist der Ansteuervektor k^p^, der die AnsteuerWerte fUr 
alle Dttsen enthSlt und folgende Bedingungen erfUllen muO: 

1. Durch k mtissen die geforderten KrSfte und Homente b. 
B ^ ^ ^ — kom 

am Raumfahrzeug ausgeUbt werden. Das heiot, die Glei- 

chung A^^^ . k^^^ „ bj^^^ erftillt werden, wobei A^ die 

Systemmatrix ist, in der die drei Kraft- und die drei ^ 

Momentenkomponenten fUr jede einzelne Diise in einer Spalte 

dieser Matrix eingetragen sind. Die Systemmatrix hat 6 Zeilen 

wegen der 6 mdglichen Freiheitsgrade in der Bewegung des 

Raumfahrzeugs, nSmlich 3 KrMfte in drei zueinander orthogonalen 

Achsen und entsprechend 3 Momente urn diese Achsen und bei 

D DUsen D Spalten. 
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2. JLqp^ « ^^opt 1 ^opt j ^opt ^^""^ """^ nichtnega. 
tive Werte enthalten, 

3. die Summe aller Ansteuerwerte muA minimal sein 



(Ziel£un)ction) , 



d.h. : 



O 

Z k 



j«1 



opt j 



Min ! 



b) Ein Simplextableau besteht aus einem Spaltenvektor der 

Basisvariablenindizes einem Zeilenvektor der Nichtbasis- 

variablenindizes % . , der Koef f izientenmatrix ^, dem 

nD -""tab 

Zeilenvektor der auf die Nichtbasisvariablen umgerechneten 

Zielf unktionskoef f izi#nten . und dem Spaltenvektor der 

■^tao *^ 

Basisvariablenwerte b^ ^: 

~tab 









• • • X ^ 










• 






• 


• 


• 


• 

^bJk 

• 




^1.1 


• 
• 


• 

aim 








... p 


... p 





Im Verfahren nach der Erfindung kommt das Simplextableau nicht 
in dieser vollen Form zum Einsatz* 

c) Filr die Umsetzung eines beliebigen 6~dimensionalen Krafte- 
Momenten-Vektors sind maximal 6 Dusen einer vorhandenen festen 
DUsenanordnung erforderlich Oder der kommandierte Krafte- 
Momenten-Vektor ist mit der vorhandenen Diisenanordnung nicht 
realisierbar. 

d) FUr jede. DUsenauswahl I von 6 aus D Dusen existiert ein 
Basisvariablenindizesvektor . Die Nummern oder Indizes 
(1,2, •••,D) der in der Auswahl I enthaltenen Diisen werden im 
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Baslsvarlablenindizesvektor Tbj^^d) zusammengefaOt und die 

Ansteuerwerte dieser DUsen sind die Werte der Basisvariablen 

des Siroplextableaus • Die Werte der Basisvariablen sind 

5 entsprechend der Reihenfolge der DUsenindizes in%^(I) im 

Basisvariablenwertevektor b^ w(^w#b ) abgelegt. Zu iedem 

""taD ""o ""Com 

BasisvariablenindizesveKtor ^j^d) gehSrt eine Basis- 
matrix Bj^(I). Die Basismatrix wird gebildet, indem man 
diejenigen Spalten der Systemmatrix A , die zu den 

10 ausgewMhlten Diisen gehdren, zu einer neuen Matrix spaltenweise 

zusammenstellt • Die Reihenfolge der Spalten in der Basismatrix 
entspricht der Reihenfolge der Indizes der DUsen im 
dazugehorigen Basisvariablenindizesvektor ^bj^d). Die 
Nichtbasismatrix B^^^j^d) wird aus denjenigen Spalten der 

15 Systemmatrix zusammengestellt , die zu den nicht in der Auswahl 

I enthaltenen Diisen gehdren. Xhre Indizes sind im Nichtbasis- 
variablenindizesvektor Ib^^j^d) zusammengefaBt . Die Reihenfolge 
der Spalten in der Nichtbasismatrix entspricht der Reihenfolge 
der Indizes der DUsen im dazugehorigen Nichtbasisvariablen- 

20 indizesvektor • Die Ansteuerwerte dieser DUsen sind Null* 

e) 1st fUr einen kommandierten KrSf te-*Momenten-Vektor b. 

""Kom 

das Optimaltableau ^ , X . ^, ^ ^ , ^ 

*^ — b.opt ' — nb.opt' —tab. opt' '*^tab.opt 

—tab opt 9®^^^**^**' dann sind die Ansteuerwerte der DUsen, die 

durch die Indizes in X ^ ^ adressiert werden, gleich Null und 

—no • opt 

25 diejenigen DUsen, die den Basisvariablen entsprechen, also 

durch die Indizes in ^ adressiert werden, enthalten als 

— b.opt 

Ansteuerwert den Wert aus der entsprechenden Zeile des opti- 

malen Basisvariablenwertevektors b^ ^ ^ ♦ Der Basisvariablen- 

— tab. opt 

indizesvektor X ^ gibt die optimale DUsenauswahl an. 
— b. opt ^ ^ 

30 f) FUr den Fall der gleichzeitigen AusUbung von alien drei 

. Kraf tkomponenten und alien drei Momentenkomponenten ist der 
Einsatz von maximal 6 der vorhandenen DUsen erf orderlich. Es 
gibt genau 



35 



Dl 



A(Cg) 



6! (D-6)l 



Mdgllchkeiten ( Kombinationen ohne Wiederholung) , 6 aus D DUsen 
auszuwahlen. Werden die Indizes in den beiden Vektoren der 
Basisvariablenindizes und Nichtbasisvariablenindizes , *X und 
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— nb' s^^^s aufsteigend sortiert, so gibt es jeweils ge- 
nau A(Cg) dieser Vektoren. FUr eine feste Dtisenanordnung gibt 
es (nur) eine bestimmte Anzahl von Kombinationen von Diisen, die 
die DUsenauswahl fUr eine Ldsung des linearen Optimierungs- 
5 problems sein konnen. Diese Anzahl ist im allgemeinen deutlich 

Kleiner als die Anzahl A(C^). 

D 

Eine wesentliche Grundlage des Verfahrens besteht darin, dafi 

stets mit sortierten Basis- und Nichtbasis variablenindizes- 

10 velctoren gearbeitet wird. Im weiteren sind die Indizes Oder 

Nummern der Diisen in den Basis- und Nichtbasisvariablenindizes- 

vektoren z.B. stets der Grofie nach sortiert* Auf diese Weise 

wird es moglich, dafi einer DUsenauswahl genau ein Indizes- 

vektor, genau eine Koef f izientenmatrix A_^ ^ « Ca^ ' , genau 

—tab tab kj ^ 

15 ein Koef f izientenvekt^r der auf die Nichtbasisvariablen 

umgerechneten Zielfunktion . und genau eine Basismatrix B. 
und Nichtbasismatrix zugeordnet werden. 

. Eine weitere wesentliche Grundlage des Verfahrens besteht 
20 darin, daB mit einem reduzierten Tableau gearbeitet wird. Das 

reduzierte Tableau besteht lediglich aus der Dilsenauswahl- 
identif ikationsnummer 1 und dem Basisvariablenwertevektor- 
FUr das Verfahren und die Arbeit mit dem reduzierten Tableau 
ist neben der Bereits tellung des Feldes der Koef f izienten- 
25 matrizen aller Opt imaltableaus und der Pivotspaltenindizes- 

matrix auch die Bereits tellung der DUsenauswahlidentif ika- 
tionsnummernmatrix, der Zeilenaus tauschindizesmatrix und der 
Basisindizesmatrix erf orderlich • 

30 Das Verfahren besteht aus zwei Bes tandteilen • Der erste, der 

Initialisierungs teil , dient zum einen der Bereit s tellung des 
Feldes der Koef f izientenmatrizen aller Optimaltableaus 

£ = [b^at kj]] = [fkjrl w<^bei [a^^ = A,,,{I) « Br\l) • B^{I) ist 
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Die DUsenauswahlidentif ikationsnummer IMuft von 1 bis 

M ( I a 1(1)M ), wobei M die Anzahl aller DUsenauswahlmoglich- 

keiten ist, die fUr einen beliebigen KrMf te-Momenten-Vektor 
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— kom vorliegende beliebige und unverSnderllche 

Diisenanordnung optimal sind. Eine Koeffizientenmatrix eines 
Optimaltableaus liegt im Sihne dieses Verfahrens dann vor, wenn 
alle Elemente im zur Koeffizientenmatrix gehdrenden Vektor der 
5 auf die Nichtbasisvariablen umgerechneten 

Zielf unktionskoef f izienten 

fi^^j^ =» ^tab'— tab' nicht negativ sind. Fiir die Elemente der 
Koeffizientenmatrix muO daher gelten: 

6 

10 1 -^^^a^ab ic j - P tab J ^ 0 ^ - '^^^'^^ 

Zum anderen werden im Initialisierungs teil des Verfahrens 
f olgende Adrefimatrizen : die Pivotspaltenindizes matrix 
JP a ^^Ilc^' die Dusenauswahlidentif ilcationsnummernma- 

15 trix U = '•^Ilc''' Zeilenaustauschindizesmatrix T = ^'^^Ijc-' 

die Basisindizesmatrix L = I^^rv^ bereitges tellt , damit auf die 
im zweiten Teil des Verfahrens benatigten Daten besonders 
schnell zugegriffen werden kann und bei der DurchfUhrung der 
Basisaustauschschritte die Arbeit mit dem reduzierten Tableau 

SO mdglich wird. Durch die DUsenauswahlidentif ikationsnum- 

mer I » ^,2, N wird jeder Dtisenauswahl , die optimal sei*n 

kann, und damit auch jeder Koef f izieatenmatrix eines Optimal- 
tableaus eine Identif ikationsnummer zugewiesen. 

15 In der Pivotspal tenindizesmatrix P = ^Pjjc^ wird fiir jede Zeile 

von jeder in £ abgelegten Koeffizientenmatrix der Index der 

dazugehorigen Pivotspalte s € {0,1,2,3, (D - 6)} abgelegt. 

Fiir eine Koeffizientenmatrix A und eine gewahlte Pivotzei- 

t a D 



30 



le 2 ist die Pivotspalte s diejenige, fiir die gilt: 
Minimum p.- n 



D.h. die zu einer Zeile einer Koeffizientenmatrix gehorige 
Pivotspalte ist diejenige Spalte, in der das Matrixelement der 
35 Koeffizientenmatrix in dieser Zeile negativ ist und der 

Quotient aus dem Element des Vektors der auf die 

Nichtbasisvariablen umgerechneten Zielf unktlonskoef f izienten in 
der Spalte und dem Betrag des entsprechenden Matrixelements aus 
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Hit Hilfe der Pivotspaltenindizesmatrlx ist flir eine Diisen-- 

auswahl I, fiir die in £ eindeutig die dazugehdrige 

5 Koef f izientenmatrix A^^j^^D abgelegt ist, und gewahlter 

Pivotzeile z die Angabe der entsprechenden Pivotspalte mit 

s(A . (I),z) » p aus der Pivotspaltenindizesmatrix moglich, 
va D X z 

ohne daO dazu die Koef f izientenmatrix wtD Oder der Vektor 

— tab 

der auf die Nichtbasisvariablen umgerechneten Zielf unktions- 
10 Roeffizienten fi^^b^"^^ benotigt wird. 

Existiert zu einer bestimmten Zeile keine Pivotspalte, so ist 

das in der Pivotspaltenindizesmatrix durch eine Null (p . » 0) 

X z 

angegeben* . 

15 

In der Diisenauswahlidentifikationsnummernmatrix U » Cu ist 
£Ur jede DUsenauswahlidentifikationsnummer I s 1,2, M und 

jede Zeile z «= 1,2,..., 6 jeder, in F abgelegten Koeff izienten- 
matrix die DUsenauswahlidentifikationsnummer derjenigen Diisen- 

20 auswahl angegeben, die nach dem Basisaustauschschritt vor- 

liegt, wenn vor dem Basisaustauschschritt die Diisenauswahl I 
vorlag und die Zeile z als Pivotzeile fiir den Basisaustausch- 
schritt gewShlt wurde. D.h., wenn zum Zeitpunkt n die Diisen- 
auswahl iCn] vorliegt und ein Basisaustauschschritt mit der 

25 Pivotzeile zCn] durchgefuhrt wird, kann die Identif ikations- 

nummer der neu entstehenden Diisenauswahl ohne Verwendung von 
A^^j^(I) mit iCn + 1] » u ^Cn] der Diisenauswahl- 

identif ikationsnummernmatrix ausgelesen werden« Dies wird 
dadurch moglich, daB durch die Kenntnis des aktuellen Tableaus 

30 (bzw. seiner Identif ikationsnummer I) und der Wahl der 

Pivotzeile die Pivotspalte implizit festgelegt ist. Zur 
Brmittlung des Folgetableaus ist also die explizlte Kenntnis 
der Pivotspalte sCn] unndtig. 

35 In der Zeilenaustauschindizesmatrix T » Ct^^l ist far jede 

— Ik 

Zeile z » 1,2, 6 von jeder der in abgelegten Koeffizien- 

tenmatrizen der Zeilenindex des Pivoteleroents im (sortlerten) 
reduzierten Tableau angegeben, das nach dem Basisaustausch- 
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schritt vorllegt, wenn vor dem Baslsaustauschschrltt die 
DUsenauswahl I vorlag und die Zeile z als Plvotzelle gewMhlt 
wurde. 

5 In der Basisindlzesmatrlx L » il^l ist fiir jede DUsenauswahl- 

Identifllcationsnummer der Vektor der Indizes der in der DUsen- 
auswahl I enthaltenen DUsen, also der zu I geh&rige Basis- 
variablenindizesvelctor ^j^d) « angegeben. Die DUsenindizes 
in den Vektoren sind der Gr6Be nach sortiert* 

0 

Durch den Initialisierungs teil des Verfahrens wird weiterhin 

eine Start-Diisenauswahlidentif ikationsnummer IC03 « I ^ ^ und 

sxart: 

die dazugehdrige invertierte Basismatrix 

'5 — b.start ^ b ^ ^start,' ^ bereitgestellt • Aus numerischen 

GrUnden 

-1 

werden alle Werte in F und B. ^ ^, die einen sehr kleinen Be- 
trag £ , (z.B*: £« lO" ) haben, zu Null gesetzt. 

10 

Im zweiten Teil des Verfahrens wird zur Realisierung der im » 
KrSf te-Momentenvektor kommandierten KrSfte und Momente die 
optimale DUsenansteuerung, d.h. die Ansteuerwerte der 
einzuschaltenden Diisen ermittelt. Nur dieser Teil des 

15 Verfahrens muB den gestellten Echtzeitanf orderungen genilgen. 

Dieser zweite Teil besteht wiederum aus zwei nacheinander 
auszuf tihrenden Teilen. Zuersi: wird ein dual zulassiges 
reduziertes Tableau als Star-ttableau ermittelt* Im zweiten Teil 
wird die erf orderliche Anzahl von Basisaustauschschritten 

10 durchgefuhrt* 

Ein dual zul^Lssiges reduziertes Starttableau wird wie folgt 

ermittelt: Fttr die Bereitstellung eines dual zulMssigen 

reduzierten Tableaus ist lediglich die Transformation des 

iS aktuellen KrSf te-Nomente-Vektors b^ auf die fttr den Start der 

Basisaustauschschritte festgelegte DUsenauswahl iCO] » I 4. . 

s t a r t 

erforderlich. Dazu dient die schon bereitgestellte invertierte 
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-1 

Basismatrix B^.g^art' ""^ ®* gilt: 
-1 

^tab. start ° °b. start * ^kom* 

5 . ■ • . 

Mit der BChon bereitges tellten Dusenauswahlidentlfikations- 

nummer fUr den Start ICOD = I ^ ^ und dem neu berechneten Vek- 

start 

'i^'tab ^9^ " ~tab start ^^'^ reduziertes und dual 

zulfisslges Starttableau CiCO] b^ ^[033 bercit. 

0 

Dann erfolgt die Ausfiihrung der Basisaustauschschrltte (zykli- 
scher Verf ahrenstell ) : 

Der folgend beschriebene Verfahrensteil muB wiederholt solange 
ausgefiihrt werden, bis die optimale Dasenauswahl mit den 
5 optimalen Ansteuerwerten gefunden wurde Oder festgestellt wird^ 

daB es fiir den aktuell vorliegenden KrSfte-Moroenten-Vektor und 
die vorliegende Diisenanordnung keine Losung gibt« 

Da von einem dual zulassigen Simplextableau in Form des redu- 
•0 zierten Starttableaus ausgegangen wird und durch die Basisaus- 

tauschschritte immer wieder nur eine Dlisenauswahl erzeugt war- 
den kann, deren Simplextableau dual zulassig ist, kann bei der 
Durchfuhrung der Basisaus tauschschritte immer wieder auf die im 
Initialisierungsteil des Verfahrens bereitges tellten Daten 
!5 insbesondere im Feld der Koef f izientenmatrizen aller optimalen 

und damit gleichzeitig aller dual zulassigen Simplextableaus 
und in der Pivotspaltenindizesmatrix zuriickgegriffen werden* 

Erster Schritt: 

10 Finde einen negativen Wert in b^ ^[n]. n in eckigen Klammern 

""tab 

bezeichnet die Nummer des Iterationsschrittes . EnthSlt b^ ^[n3 
keinen negativen Wert, liegt die LSsung des linearen Optimie- 
rungsproblems vor, und es kann der Vektor der Dttsenan- 
steuerwerte jL^p^ aufgebaut werden. EnthMlt b^^b'"'^-' "^ST^^i^^ 
5 Werte, dann wird eine dieser Zeilen als Pivotzeile zCn] 

gewMhlt. Es ist gilnstig, die Zeile des betragsgrSBten negativen 
Elements in k^^j^Cn] als Pivotzeile zu wfihlen. 
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Zweiter Schrltt: 

Mit Kenntnis der DUsenauswahlidentifikationsnummer I[n], und 
der im ersten Schrltt gewahlten Pivotzeile zCn] wird aus der 

5 Pivotspaltenlndizesmatrlx JP der dazugehdrige Pivotspaltenindex 

mit sCn] « p T-r-^i-r n ausgelesen. 1st in P kein Pivotspal- 
tenindex angegeben (p jf^lzCn] ° 0),muU eine andere Pivotzeile 
ausgewShXt -werden. 1st zu keiner der mdglichen Pivotzeilen in 
der Pivotspaltenindizesmatrix eine Pivotspalte indiziert, 

0 existiert keine L6sung des linearen Optimierungsproblems « Sind 

Pivotzeile und Pivotspalte ermittelt, wird zum dritten Schritt 
ilbergegangen. 

Im dritten und im vierten Schritt wird der eigentliche Basis- 
IS austauschschritt durch Ermittlung des nachsten reduzierten 

Simplextableaus ausgefUhrt: 
Dritter Schritt: 

Mit Hilfe der Koef f izienten in der Pivotspalte sCn] von wCn] 

^tab 

wird aus b . Cn] der neue Basisvariablenwertevektor b^ . 

"^tab r-tab«neu 

^0 berechnet. Die benStigten Elemente des Pivotspaltenvektors a 

aus A^ajjCj^^ konnen direkt aus dem Feld der Koef f izienten- 

matrizen aller Optimaltableaus mit a^ ^ , r iCn] « £, r i^r t 

tab k sCn] k sLnjXLnJ 

ausgelesen werden. 

Die Elemente des neuen Basisvariablenwertevektors werden wie 
IS folgt berechnet: 

Der neue Wert in der Pivotzeile zCn] wird berechnet nach: 

b - ^1^1 
tabjteu zln] T , 

JO wobei ^^1^2 sCn] iCn] entsprechend indizierte Pivotelement 

aus A^ab^^^ ist • 

Die neuen Werte in den anderen Zeilen werden berechnet nach: 

^ " ^rab k C«] - ^a6^ z[n] ' fks[n]I[nl 
wobei: {k = 1(1)6) aber [k 9^ z[n]) 

Da der neuberechnete Basisvariablenwertevektor b^ . zum 

*~taD*neu 
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flktlv angenommenen neuen sortierten Basisvariablen- 
indizesvektor 'X^Cn +1] passen muB, ist Im allgemeinen der 
neuberechnete Baslsvarlablenwertevektor entsprechend 
umzuordnen. Dazu sind Im neuberechneten Baslsvarlablen- 
wertevektor b. ^. _ ledlgllch die Werte zwlschen der aktuellen 
— tab • neu ^ 

Plvotzelle z[n] und der mlt Hllfe von I[n] und zCn] aus der 

Zellenaustauschlndlzesmatrlx T = [t^ ] auszulesenden Austausch- 

zelle Zj.^,, A^uf ^ -.r 1 r T zyklisch zu tauschen« Die 
taudcn ILnJ zLnj ^ 

Austauschzelle z^ ^[n] 1st die Zelle, wohln der Wert des 

tausch ' 

neuberechneten Basisvarlablenwertevektors plaziert wlrd, der in 

b. in der Plvotzelle zCn] steht. Nach dlesem Umordnen des 
**tao«neu 

neuberechneten Basisvarlablenwertevektors llegt der fUr den 
f olgenden Baslsaustauschschritt ben&tlgte Baslsvarlablen- 
wertevektor ^^g^jjCn + 1] vor. 

Vierter Schritt: 

Mlt Kenntnis von iCn] und der Plvotzelle zCn] wird mlt 

iCn + 1 3 » u^r_^ - . die Dusenauswahlidentlflkatlonsnummer 
XLnj ZLnj 

ermlttelt, die zum neuberechneten und nach auf steigenden 
Basisvariablenindizes sortierten Baslsvarlablenwertevek- 
tor ^^gjjCn + 13 gehttrt. Der Baslsaustauschschritt ist damit . 
vollzogen, das neue reduzierte Tableau [I[n +13 b^^^^^Cn + 133 
llegt vor, und es kann die Losung des linearen Optimlerungs- 
problems mlt dem ersten Schritt des zykllschen Verf ahrensteils 
fortgesetzt werden* 

Bricht das Verfahren mlt elnem Baslsvarlablenwertevektor, der 

keine negativen Werte mehr enth£llt, ab, so llegt der optlmale 

6-dimensionale Baslsvarlablenwertevektor b^ . = 

—tab .opt —tab end 

vor, und es kann mlt Hllfe von iCn .3 aus der Basislndlzes- 

eno 

matrix L der zum aktuellen Baslsvarlablenwertevektor gehSrende 

Vektor der Basisvariablenindizes mlt VbCn .3 = ^xr^ i 

end iLn j 

end 

ausgelesen werden. Hit Kenntnis der beiden Vektoren b. .[n .3 

""tab ena 

— b'^'^end*' ^^^^ optimale D-dim^nsionale Ansteuervektor 

fiir die DUsen li^p^ aufgebaut. 

FUr die LSsung des linearen Optimlerungsproblems fUr 
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unterschledliche KrSf te*-Momenten-Vektoren sind im allge- 

meinen auch unterschledlich vlele Baslsaustauschschrlt* 

te n^^, € {0, 1, n ) erforderllch. Bei Simulatlonsrech- 

end max 

nungen mit elner Diisenanordnung von 12 Diisen wurden in sehr 

seltenen FMllen bis zu n » 10 Basisaustauschschritte durch- 

max 

geftihrt, obwohl im Durchschnitt nur ca. 3,5 Basisaustausch- 
schritte notig waren. Bei geeigneter Wahl des Starttableaus 
^^start — tab.start"' Jedoch nur in ca. 0,5 96 aller Fillle 

mehr als 5 Austauschschritte notwendig. Im allgemeinen bleibt 
es ohne merlcliche Konsequenzen fiir die Anwendung, wenn die 
Anzahl der maximal auszufUhrendeh Austauschschritte auf einen 
Wert begrenzt wird, bei dem nahezu alle FSlle vollstandig 
bearbeitet werden und in den Fallen, in denen doch noch 
negative Werte im Basisvariablenwertevektor stehen, diese zu 
Null gesetzt werden. Auf diese Weise kann der Worst-Case- 
Rechenaufwand weiter gesenkt werden. Anwendungen dieser 
Verfahrensmdglichkeit bei der Simulation in geschlossenen 
Regelkreisen zeigten keine bemerkbare Verschlechterung der 
Regelgiite Oder eine Erh6hung des Treibstof f verbrauchs . 

Enthalt der Basisvariablenwertevektor nach der Ausfiihrung der 
maximal erlaubten Anzahl von Basisautauschschritten noch 
negative Werte, ist es prinzipiell auch moglich, daB keine 
LSsung des linearen Optimierungsproblems existiert. Dieser Fall 
kann dann au£treten, wenn eine Dusenanordnung vorliegt, die 
nicht beliebige Krfifte und Momente aus dem Bereich aller 
zugelassenenen Kraf te-Momenten-Vektoren erzeugen kann. Damit in 
diesen FMllen keine zu groBen Fehlkrafte und Fehlmomente 
erzeugt werden, kann in Abhangigkeit vom Betrag der negtiven 
Elemente im Basisvariablenwertevektor entschieden werden, ob 
die negativen Elemente zu Null gesetzt werden und der Ansteuer- 
vektor fiir die Dusen aufgebaut wird Oder generell keine An- 
steuerung der DUsen erfolgt. 

kommen bei der Abarbeitung des Simplexalgorithmus entartete 
Optimaltableaus vor, dann ist es theoretisch m5glich, daB der 
Algorithmus in einen Endloszyklus lauft und nicht konvergiert. 
In Simulationsrechnungen mit unterschiedlichen DUsenanordnungen 



wo 98/49058 



PCT/EP98/0159S 



1st dieser Fall allerdings noch nle elngetreten und wlrd 
deshalb auch nicht expllzit berUcksichtlgt • FUr den dualen 
Slmplexalgorithmus kann man durch die gewShl-te Dtisenanordnung 
sicherstellen, daB keine dual entarteten Optimaltableaus 
existleren. FUr den primalen Slmplexalgorithmus, also das 
Auftreten primal entarteter Optimaltableaus, ist das so nicht 
mSglich, 

Soil bei der Ermittlung der treibs tof f minimalen Triebwerks- 
ansteuerung der unterschiedliche spezifische Treibstoff- 
verbrauch , j«l(1)D der einzelnen Triebwerke beriicksichtigt 
werden, liegt die Zielfunktion 

D 

Z c k « Min ! 

vor. ^ 

Die Berttcksichtigung der Gewichtskoef f izienten ist in dem 
Verfahren filr den on-line auszufUhrenden Verfahrensteil ohne 
Mehraufwand moglich. 

F£Lllt eine Diise aus , andert sich die DUsenanordnung im Sinne 
dieses Verfahrens. Das Verfahren auf die neue DUsenanordnung 
anzuwenden eirfordert, daB alle Daten, die im Initiali- 
sierungsteil des Verfahrens bereitgestellt wurden, fiir die neue 
DUsenanordnung bereitgestellt werden. 

Da sich die neue DUsenanordnung nur durch den Wegfall einer 
DUse von der alten DUsenanordnung unterscheidet , gibt es 
M5glichkeiten, die neuen Daten durch eine Aktualisierung der 
alten Daten zu ermitteln* Durch die Bereitstellung eines 
erweiterten Feldes von Koeffizientenmatrizen, in dem auch 
diejenigen Koeffizientenmatrizen abgelegt sind, die durch den 
Ausfall von DUsen Koeffizientenmatrizen von Optimaltableaus 
werden, und entsprechend erweiterter AdreBmatrizen, kann der 
Aktualisierungsauf wand stark reduziert werden und fUr bestimmte 
Matrizen ganz entfallen. 

SchlieBlich besteht die Mdglichkeit, die Daten fttr alle 
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DUsenanordnungen, die durch den Ausfall einer DUse zu erwarten 
sind, vorab bereitzustellen , und bei Ausfall einer bestimmten 
DUse dann auf die entsprechenden neuen Daten umzuschalten. Bis 
zum Ausfall einer nfichsten Diise miissen fur die dann zu 
ervrartenden DUsenanordnungen die Daten wiederum bereitgestellt 
werden, urn dann durch Umschalten auf den nSLchsten DUsenausfall 
reagieren zu Icdnnen. 

Auch bei einer Verschiebung des Schwerpunktes des Raumfahr- 
zeuges kommt es zu einer VerMnderung der Diisenanordnung im 
Sinne dieses Verfahrens* Iro allgemeinen fiihrt das zu einer 
VerMnderung der Systemmalirix* Daher miissen alle von ihr abge- 
leiteten Daten neu bereitgestellt werden. 

Das zuvor beschriebene Verfahren kann zur weiteren Einsparung 
von Rechenauf wand und Speicherplatzbedarf noch modifiziert 
werden: 

- In vorans tehender Beschreibung wurde stets von der Steuerung 
aller 6 Freiheitsgrade des Raumf ahrzeugs ausgegangen* Das 
Verfahren ist auch auf die Steuerung von 1 bis 5 Freiheits- * 
graden anwendbar. Es ist lediglich die Streichung bzw* die 
Nichtberucksichtigung der den nicht zu steuernden Freiheits- 
graden entsprechenden Zeilen in der Systemm'atrix und der von 
ihr abgeleiteten Vektoren und Matrizen erf orderlich, 
Desgleichen ist die Steuerung einer sich standig Mndernden 
Teilmenge aller 6 Freiheitsgrade moglich. Liegt eine schiefe 
Diisenanordnung vor, bleiben jedoch die nicht gesteuerten 
Freiheitsgrade im allgemeinen nicht unbeeinf luBt . 

- Miissen am Raumfahrzeug nicht beliebige KrAfte und Momente 
ausgeiibt werden, ist der Einsatz einer reduzierten Diisen- 
anordnung moglich. Das vorliegende Verfahren ist auch fiir 
derartige Anordnungen anwendbar. Es sollte nur gewSLhrleistet 
sein, dal) die kommandierten Kr^ifte und Momente des KommandQs 
durch die vorliegende Diisenanordnung realisierbar sind. Ist das 
Kommando nicht realisierbar, gibt das Verfahren an, daB keine 
Ldsung des linearen Optimierungsproblems existiert. 
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- In Abhangigkeit von der vorliegenden Dusenanordnung kommt es 
vor, daB von den im Feld der Koef f izientenmatrizen aller 
Optimaltableaus abgelegten Koef f izientenmatrizen bestimmte 

5 Spalten nie als Pivotspalte ausgewMhlt werden. Nach Berechnung 

der Pivotspaltenindizesmatrix P kdnnen die nicht benStigten 
Spalten aller Koef f izientenmatrizen einfach ermittelt werden. 
Es ist moglich, ein Feld aller relvanten Pivotspaltenvektoren V 
zu erstellen, in dem nur die Spalten der Koef f izientenmatrizen 

10 aller Optimaltableaus abgelegt sind, die tatsachlich als Pivot- 

spaltenvektoren benotigt werden. Dadurch kann eine signifikante 
Speicherplatzersparnis erreicht werden* Die Adressierung der 
Matrixspalten im Feld aller relevanten Pivotspaltenvektoren 
erfordert keinen erhohten Rechenauf wand. Jede Spalte hat im 

15 Feld aller relevan-ten Pivotspaltenvektoren ihren Spaltenindex, 

Wird in einer verMnderten Pivotspaltenindizesmatrix, der 
Pivotspalteniden-tif ikationsnummernmatrix 0> = C^tv^ nicht der 
Spaltenindex eines Pivotspaltenvektors in einer 
Koef f izientenmatrix angegeben, sondern gleich die 

20 Pivotspaltenidentif ikationsnummer als Spaltenindex r im Feld 

aller relevanten Pivotspaltenvektoren V, dann kann im zweiten 
Schritt des zyklischen Verf ahrensteils mit r » o^j^ gleich der 
Spaltenindex im Feld aller relevanten Pivot-spaltenvektoren V « 
CVj^j.] angegeben werden. Die Ermittlung der Elemente des neuen 

25 Basisvariablenwertevektors im dritten Schritt des zyklischen 

Verf ahrensteils andert sich dann wie folgt: 

. _ ^tab zfw] t^I 

wobei ^2[n] r[n] entsprechend indizierte Pivotelement aus 

A^g^jjjCn] ist* Die neuen Werte in den anderen Zeilen werden 
berechnet nach: 

wobei: (k = 1(1)6) aber {k ^ z[n]) 
Die weiteren Verf ahrensbestandteile bleiben unverandert* 

- Liegt ein Feld aller relvanten Pivotspaltenvektoren V vor. 



30 



35 



BNSOOCiD: <WO ^9849058A1J_> 



19 



wo 98/49058 



PCT/EP98/01595 



kann insbesondere bei DUsenanordnungen, die bestimmte 
Symmetrlen aufweisen, festgestellt werden, daB relevante 
Pivotspaltenvektoren unterschiedlicher Optimaltableaus 
numerisch identisch sind. Da das mehrmaiige Ablegen gleicher 
Pivotspaltenvektoren nicht notwendig ist, kann erneut weiterer 
Speicherplatz eingespart werden. Beim sukzessiven Aufbau des 
Peldes aller relevanten Pivotspaltenvektoren im Initiali- 
sierungsteil des Verfahrens kann uberpriift werden, ob ein 
numerisch identischer Pivot spaltenvektor zum gerade ermittelten 
Pivotspaltenvektor schon in V abgelegt wurde. 1st dies der 
Fall, kann das erneute Ablegen dieses relevanten 
Pivotspaltenvektors entfallen und in der zu V gehorigen 
Pivotspaltenidentifikationsnummernmatrix 0. wird die 
Pivotspaltenindentifikationsnummer ( Spaltenindex) des schon im 
Feld aller relevanten Pivotspaltenvektoren abgelegten 
Pivotspaltenvektors eihgetragen. 

~ Besteht das Erfordernis^ den Speicherbedarf des Verfahrens 
noch vreiter zu reduzieren, kann eine Verfahrensvariante 
eingesetzt werden, bei der es nicht erforderlich ist, alle 
relevanten Pivotspaltenvektoren, sondern nur noch ReprHsen- ^ 
tantenvektoren im Rechnerspeicher abzulegen. ReprSsentanten- 
vektoren sind relevante Pivotspaltenvektoren, aus denen mit nur 
minimalem Rechenauf wand weitere relevante Pivotspaltenvektoren 
berechnet werden konnen. 

T 

Jeder Pivotspaltenvektor a^ = Ca^a^ sl^ a^] gehort zu 

einer bestimmten Basis, z«B. 

^wdi ) « -I "'^i o ••• ''^ ^ 3*^ und zu einer 

b 1 k«il.2 s»i p.n 

Nichtbasisvariablen, z.B.X ,^ '^• • Wenn a. verschieden von 

nb r X n* 

Mull ist, lant sich aus a(I^,>.^) « ^^1^2 *i ***^n^ 

derjenige Pivotspaltenvektor 

^neu " •i<^2'^s> " ^ *neu. 1 *neu. 2 ^neu.i ^neu.n^ 



berechnen, der zur Basis 
Jo 2 k.11.2 r.i p.n 
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Nlchtbasisvariablen X o X gehdrt. So kann aus a(I,, ) mit 

no s ^ I r 

der Angabe des Austauschlndexes i der neue Pivotspaltenvektor 
aCl^, berechnet werden. 

5 Auf diese Weise konnen aus einem z.B. sechsdimensionalen 

Pivotspaltenvektor bis zu sechs weitere Pivotspaltenvektoren 
berechnet werden. FUr die Anwendung im Verfahren zur treib- 
stoffminimalen Diisenans teuerung kann man sich jetzt auf das 
Speichern derjenigen relevanten Pivotspaltenvektoren als 
10 ReprMsentantenvektoren beschrSnken, die erforderlich sind, urn 

alle relevanten Pivotspaltenvektoren wie angegeben zu 
berechnen. 

Piir die Ermittlung einer mSglichst kleinen Anzahl von 
15 ReprSsentantenvektoren lassen sich sehr gilnstig Genetische 

Algorithmen als Optiomierungsverfahren einsetzen. 

Fiir eine geeignete Adressierung der ReprSsentantenvektoren in 
der Matrix aller Reprasentanten sind statt der 

20 Pivotspaltenindizes - bzw. 

Pivotspaltenindentifikationsnummernmatrix einmal eine 
ReprMsentantenidentifikationsnummernmatrix R und weiterhin eine 
Austauschindizesmatrix erforderlich. In der 
Reprasentantenidentif ikationsnummernmatrix R wird ftir jede 

25 DUsenauswahlidentif ikationsnummer und jede Pivotzeile die 

Identif ikationsnummer des Reprasentantenvektors angegeben, der 
fiir die Berechnung des tatsachlich benotigten 
Pivotspaltenvektors verwendet werden kann* In der 
Austauschindizesmatrix S wird fur jede 

30 Dusenauswahlidentif Ikationsnummer und jede Pivotzeile der 

Zeilenindex des Aus tauschelements angegeben, um aus dem 
ReprSsentantenvektor den tatsachlich bendtigten Pivotspalten- 
vektor berechnen zu. kdnnen. Steht an einer Stelle in der 
Austauschindexmatrix eine Null, bedeutet dies, daB der 

35 adressierte ReprMsentantenvektor unmittelbar der erf orderliche 

Pivotspaltenvektor ist und daher die Umrechnung in den aktuell 
bendtigten Pivotspaltenvektor entfMllt. 
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Da stets mit nach Basisvariablenindlzes sortierten 
Pivotspaltenvektoren gearbeitet wird, muB der zun&chst aus dem 
ReprSsentantenvektor berechnete PivotspaltenveJctor noch 
umsortiert werden (z.B. nach aufsteigenden 
5 Basisvariablenindlzes), urn benutzt werden zu k&nnen. Wlrd die. 

Sortierung der GroOe der Baslsvarlablen vorgenommen, glbt es 
bei elnem sechselementlgen Vektor genau 30 M6gllchkelten der 
Umsortlerung, urn die erf or-derllche Ordnung wieder 
herzustellen* Bei einer erforderlichen Umsortierung wlrd stets 
10 ein Vektorelement an einer Position im Vektor entnommen und auf 

eine andere Position geschrieben* Die Vektorelemente zwischen 
den beiden Positionen werden entspre-chend verschoben. 

In einer Sortlerldentif Ikationsnummernmatrlx S Cfl j wlrd 
15 £Ur jeden ReprSsentanten und Jede Austauschzelle die Sortler- 

identifikationsnuniiner angegeben, nach der der Vektor umgeordnet 
werden muB, um den benStlgten Plvotspaltenvektor zu erhalten. 

Aus elneiii ReprMsentantenvektor k5nnen durch Baslsaustausch 
20 und/oder oben beschriebene Umordnung bis zu 7 x 31 » 217 

Vektoren geblldet werden. Werden alle dlese Mdgllchkelten 
berilckslchtlgt , kann es aufgrund von Symmetrien In der 
Triebwerksanordnung vorkommen, daB wenlger Reprasentanten- 
vektoren erforderllch sind, als wiirden nur die sieben aus einem 
25 ReprMsentantenvektor gebildeten relevanten Pivotspaltenvektoren 

in Betracht gezogen. 

Eine weitere Ver£ahrens variante ergibt sich daraus, daB auf den 
Baslsaustausch am ReprMsentantenvektor verzichtet wlrd und nur 

30 die 30 Umordnungsmogllchkelten des Reprasentantenvektors 

berlickslchtlgt werden. Auf dlese Welse reprMsentlert ein 
ReprMsentantenvektor 31 Vektoren. Aufgrund von Symmetrien In 
der Triebwerksanordnung k6nnen auch auf dlese Helse mehrere 
relevante Pivotspaltenvektoren generlert werden. Im Verglelch 

35 zur letzten Verf ahrensvarlante werden im allgemelnen jedoch 

mehr Reprasentantenvektoren benotlgt, jedoch entfMllt 
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glelchzeitig der Rechenaufwand fiir den Basisaustausch am 
Reprfisentantenvektor, und es ist die Austauschlndizesmatrix 
nicht erf order lich. 
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Patentansprilche 

1 • Verfahren zur Insgesamt trelbstoffminimalen rechnerge- 

stiitzten Ansteuerung einer bellebig festgelegten Anzahl und 

beliebig festgelegten Anordnung von DUsen an eineni 

Raumf ahrzeug, 

dadurch gekennzelchnet, 

- da0 ein treibstof f minlmaler Ansteuervektor (k ^) mit 

opx 

ausschlieBlich nicht negativeh Ansteuerwerten fiir alle 
vorgegebenen DUsen mit Hllfe des dualen Slmplexalgorithmus 
gebildet wird, 

- da& dazu jeder DUsenauswahl aus der Gesaml:zahl von Diisen am 

Raumfahrzeug, die die optimale DUsenauswahl fUr die treib- 

stoffminimaie Realisierung eines beliebigen Kommando- 

vektors (b, ) darstellt. eine DUsenauswahliden- 
— kom ^ ' 

tif ikationsnummer (1) zugeordnet wird, 

* daB fUr die Ermittlung der optimalen Ansteuerung zunMchst 

der Kommandovektor (b, ) der die auszuUbende Kraft 

— kom 

und/oder Momentenkomponenten enthfilt, durch Multiplikation 

mit einer Transformationsmatrix ((B ^ ^)^^) in einen 

start 

Startvektor (b=b ^ ^) transformiert wird, 

—start 

daB danach der duale Simplexalgorithmus durchgefUhrt wird, 
bis der Vektor der Ansteuerwerte (b) nur nichtnegative 
Ansteuerwerte enthMlt oder aber die Nichtexistenz einer 
L5sung angezeigt wird, 

- und daB dabei der duale Simplexalgorithmus realisiert wird, 

indem vorab alle Daten, die fUr die Abarbeitung des dualen 
Simplexalgorithmus notwendig sind und sich off-line (vor- 
ab) berechnen lassen, ermittelt werden, in Datenfeldern 
im Rechnerspeicher abgelegt werden und beim Einsatz des 
dualen Simplexalgorithmus wahrend des on-line-Betriebes 
durch geeignete Adressierungen aus dem Rechnerspeicher 
abgerufen und in der Rechnung verwendet werden, wobei nur 
mit einem reduzierten Simplextableau gearbeitet wird, in- 
dem ledigl.ich die DUsenauswahlidentif ikationsnummer und 
der dazugehSrige Vektor der Ansteuerwerte (b) gefUhrt 
werden. 
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2. Verfahren nach Anspruch 1, 
dadurch gekennzelchaet, 

dafi die Transformatlonsmatrlx ((B ^ ^)^^) die invertierte 

sxarr 

Baslsmatrix derjenlgen Diisen darstellt, die zur Startdttsen- 

5 auswahl gehdren, diese DUsenauswahl fttr einen beliebigen 

Kommandovektor optimal ist, diese DUsenauswahl eine Start- 

DUsenauswahlidentif ikationsnummer (I . . ) zugeordnet ist und 

s ^a irc 

der Startvektor (b«b^^^^^) die Ansteuerwerte dieser Dusen 
enthfilt. 

10 

3. Verfahren nach einem der Anspriiche 1 Oder 2, 
dadurch gekennzeichnet , 

- dafi jede Diise am Raumfahrzeug eine Identif ikationsnummer 

erhMlt, und die Werte im Vektor der Ansteuerwerte fiir eine 
15 bestimmte DUsenauswahl stets in einer festgelegten Ordnung 

anhand der Identifikationsnummern der DUsen sortiert 
werden, 

- dafibei der Abarbeitung des dualen Simplexalgorithmus nach 

der Ermittlung einer Pivotzeilennummer und mit Kennthis 
20 der DUsenauswahlidentifikationsnummer der aktuell vor- 

liegenden DUsenauswahl die entsprechende Pivotspalten- • 
nummer aus einer vorab bereitgestellten Pivotspaltenin- 
dizesroatrix ausgelesen wird, 

- daB das Nichtvorhandensein einer Pivotspal tennummer zu der 
25 vorliegenden DUsenauswahlidenti£ikationsnummer und zu 

der gewahlten Pivotzeile durch einen besonderen Eintrag 
in der Pivotspaltenindizesmatrix eindeutig angezeigt wird, 
-* und daB die Nichtexis tenz einer Losung dadurch angezeigt 
wird, daB zu keiner moglichen Pivotzeile in der Pivotzei- 
30 lenindizesmatrix eine Pivotspaltennummer angegeben ist. 

4. Verfahren nach einem der AnsprUche 1 bis 3, 
dadurch gekennselchnet, 

- daB der Pivotspaltenvektor der Koef f izientenmatrix des 

35 aktuellen Simplextableaus , dessen Werte fttr die Berechnung 

neuer Ansteuerwerte fUr die DUsen einer neuen DUsenauswahl 
benotigt werden, mit Kenntnis der DUsenauswahl vor dem 
Basisaustauschschritt und der gewMhlten Pivotspaltennummer 
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au8 elnem vorab bereitgestellten Feld ailer moglichen 
optimalen Koeff izientenmatrlzen ausgelesen wird, 

- daB der Vektor der Ansteuerwerte ftir die Diisen der Diisen- 

auswahl, die nach dem Basisaustauschschritt vorliegt, 
berechnet wird und die Werte in diesem Vektor entsprechend 
der Zuordnung zu den DUsen und entsprechend der festgeleg- 
ten Ordnung der Identifikationsnummern fUr die Dttsen 
sortiert werden. 

5* Verfahren nach einem der Anspriiche 1 bis 4, 
dadurch gekennzeichnet , 

- daa mit Kenntnis der Dusenauswahlidenti£ikation8nummer vor 

dem Basisaustauschschritt und der gewShlten Pivotzeilen- 
nummer die Dusenauswahlidentif ikationsnummer der jenigen 
DUsenauswahl , die nach dem Basisaustauschschritt 
vorliegt , aus einer DUsenauswahlidentifikationsnummern- 
matrix ausgelesen wird 

- und dafi nach Voriiegen eines Vektors der Ansteuerwerte, der 

nur noch nicht negative Werte enthMlt, mit Kenntnis der 

dazugeh5rigen DUsenauswahXidentif ikationsnummer der Vektor 

der Dtisenindizes aus einer Basisvariablenindizesmatrix * 

ausgelesen wird und die mit ihrer DUsenidentifikations- 

nummer indizierten DUsen den entsprechenden Wert aus dem 

Vektor der Ansteuerwerte (b) erhalteh sowie die iibrigen 

DUsen den Ansteuerwert Null zugewiesen bekommen und damit 

der Ansteuervektor (k ^) bestimmt ist. 

"~opt 

6. Verfahren nach einem der Anspriiche 1 bis 5, 
dadurch gekennzeichnet , 

daB nur diejenigen Spalten von Koeff izientenmatrizen dual 
zulassiger Simplextableaus in einer Matrix aller relevanten 
Pivotspaltenvektoren abgelegt werden, die bei der Abarbeitung 
des Verfahrens tats^lchlich als Pivotspaltenvektoren ausgewahlt 
warden kdnnen/ 

und daB identische relevante Pivotspaltenvektoren nur einmal in 
dieser Matrix abgelegt werden. 
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7. Verfahren nach Anspruch 6, 
dadurch gekennseichnet, 

- daO £Ur die Adresslerung der relevanten Pivotspalten- 
5 vektoren jedem Plvotspaltenvektor in der Matrix aller 

relevanten Pivotspaltenvektoren eine Pivotspalteniden- 
tifikationsnummer zugewiesen wird und diese in eiher 
Pivotspaltenidentif ikationsnummernmatrix derart eingetra*. 
gen werden, daO bei Kenntnis der aktuellen DUsenauswahl- 
0 identif ikationsnummer und der aktuell gewMhlten Pivotzei- 

lennummer die Pivotspaltenidentifikationsnummer des aktu- 
ell ben6tigten Pivotspaltenvektors auslesbar ist, 

- daB das Nichtvorhandensein einer Pivotspalte zu einer 

gewMhlten Piv.otzeile durch einen besonderen Eintrag in der 
5 Pivotspaltenidentifikationsnummernmatrix eindeutig 

angezeigt wird, 

- und daB die Nichtexistenz einer L5sung dadurch bestimmt 

wird, daB zu keiner moglichen Pivotzeile in der Pivot- 
spaltenidentifikationsnummernmatrix eine Pivotspalten- 
0 identifikationsnummer angegeben ist. 

8. Verfahren nach einem der Anspriiche 3 bis 7, 
dadurch gekennzelchnet , 

daB statt der relevanten Pivotspaltenvektoren nur bestimmte 
:5 Pivotspaltenvektoren als Reprasentantenvektoren im 

Rechnerspeicher abgelegt werden, wobei die 
Reprasentantenvektoren eine mdglichst kleine Anzahl von 
Pivotspaltenvektoren darstellen, aus denen die tatsachlich 
ben5tigten relevanten Pivotspaltenvektoren ermittelt werden, 
;0 und daB jeder Reprasentantenvektor eine 

Reprasentantenidentifikationsnummer zugewiesen bekommt und zur 
Adressierung der ReprMsentantenvektoren eine 

ReprMsentantenidentif ikationsnummernmatrix verwendet wird, in 
der £Ur jede Dtisenauswahlidentifikationsnummer und jede 
;5 Pivotzeilennummer die ReprMsentantenidentifikationsnummer 

desjenigen ReprMsentantenvektors eingetragen ist, der filr die 
Berechnung des erforderlichen Pivotspaltenvektors verwendet 
wird. 
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9« Verfahren nach Anspruch 8, 
dadurch gakennselclinet , 

daG In elner Austauschmatrix fUr jede DUsenauswahl- 
identlflkationsnummer und jede. Pivotzellennummer der 
Zeilenlndex der auszutauschenden Baslsvariablen als 
Austauschelement angegeben ist, um aus dem ReprSsentantenvektor 
den benS.tigten Pivotspaltenvektor zu berechnen, 

10. Verfahren nach einem der Anspriiche 1 bis 9, 
dadurch gekennzeichnet , 

daB £Ur die M&glichkeit der Veranderung der Dtisenanordnung 
durch zu erwartende Diisenausfalle die entsprechenden Daten 
vorab berechnet und gespeichert werden und nach Feststellen 
eines D(isenaus£alls automatisch durch den Rechner au£ Verwen- 
dung der entsprechenden neuen Daten umgeschaltet wird. 

11. Verfahren nach Anspruch 10, 
dadurch gekennzelchnet , 

daB die nach einem Diisenausf all erf orderlichen neuen Daten 
durch eine entsprechende Korrektur und Aktualisierung der bis 
zum Diisenausfall verwendeten Daten bereitgestellt werden. 

12. Verfahren nach einem der Anspriiche 1 bis 11, 
dadurch gekennseiehnet , 

daB das Verfahren unabhangig von der Dimension des Kommando<* 

vektors (b, ) angewendet wird und dabei fiir die Berechnung der 
"~Jcom 

vorab bereitgestellten Daten die nichtauszuUbenden Komponenten 
aus der vollstandigen Systemmatrix der Diisenanordnung 
eliminiert werden. 
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wahrand der Internatlonaten Recherche konaultlefte elektroniache Datenbank (Name der Datenbank und evtl. verwendete SuchbegrWe) 



C. ALS WESENTLICH AN6ESEHENE UNTERLAGEN 



Kategorie^ Bezeichnung der Verdffentllchung. sowelt erfordarlk:h unter Angabe der in Betracht kommandan Tetta 



Betr. Anspruch Nr. 



WIKTOR P J: "MINIMUM CONTROL AUTHORITY 
PLOT: A TOOL FOR DESIGNING THRUSTER 

SYSTEMS" 

JOURNAL OF GUIDANCE AND CONTROL AND 
DYNAMICS, 

Bd. 17, Nr. 5, September 1994, 

Selten 998-1006, XP000493390 

1n der Anmeldung erwahnt 

siehe Selte 1001, Spalte 2, Absatz 8 - 

Selte 1002, Spalte 1, Absatz 2 



1.2 



US 5 130 931 A (PALUSZEK MICHAEL A 
14. Jul i 1992 

siehe das ganze Dokument 



ET AL) 



1.2 



-/-- 



Weitere Ver6ffentlichungen sind der Fortsetzung von Feld C zu 
entnehmen 



Siehe Anhang Patentfamilie 



Besondere Kategorien von angegebenen Ver6ffentllchungen 
"A" Vardffenttlehung, die den allgemeinen Stand der Technik detlnleit. 
aber nicht ala beaondera bedeutsanrtanzusahen 1st 

"E" Sltares Dokument. daa jedoch erst am oder nach dam IntomattorMilen 
Anmaldedatum verdffentlicht worden Ist 

"L" Verdffentllchung. die geelgnet let, einen Prioritdtsanspruch zweifelhaft er- 
schelnen zu laasen. Oder durch die das Verdffentlichungsdatum einer 
anderen im Recherchenbericht genannten Verdffentlichung belegt werden 
soti Oder die aus einem anderen beaonderen Qrund angegeben let (wto 
ausgofOhrt) 

"O" Veroftentiichung, die sich auf eine mundltche Offenbamng. 

eine Benutzung. eine Ausstellung Oder andare MaOnahmen bezieht 
'P" Verdffentlichung, die vor dam tntemationaten Anmekiedatum. aber nach 

dam beanspruchten Prioritfltsdatum verdffentlicht worden let 



T* SpStere Verdffentlichung. die nach dam Internationalen Anmaldedatum 
Oder dem Priorit^sdatum verdffentircht worden ist und mit der 
Anmeklung nicht koilUtert. sondarn nur zum VerstSndnts des der 
Erflndung zugrundellegenden Prinzlps Oder der ihr zugrundellegenden 
Theorie angegeben ist 

"X" Verdffantltehung von besonderer Bedeutung; die beanspruchte Erfindung 
kann allein aufgmnd dieser Verdffentlichung nicht als neu oder auf 
erfinderischer Tdtigkeit beriihend betrachtet werden 

"Y" Verdffentltehung von besonderer Bedeutung; die beanspruchte Erfindung 
kann nicht ais auf erfinderischer Tatlgkeit beaihend betrachtet 
werden. wenn die Verdffentllchung mit einer oder mehreren anderen 
Verdffentllchungen dieser Kategorie in Verbindung gebracht wlrd und 
diese Verbindung fur einen Facnmann naheiiegend Ist 

Verdffentlichung, die mtgUod deraelbenPatentfamiile ist 



Datum des Abachlusses der internationalen Recherche 



3.JU11 1998 



Absendedatum des Internationalen Recherchenberfchta 



16/07/1998 



Name und Postanschrift der Internationalen Recherchenbehdrde 
Europdisches Patentamt« P.B. 5816 Patentlaan2 
NL - 2280 HV Rijswijk 
Tel. (4^1-70) 340-2040, Tx. 31 651 epo nl. 
Fax: (431-70) 340^16 



BevollmAchtigter Bediensteter 

Calvo de No, R 



FbrnlUatt PCT/ISA^O (Batt 2) (JuB 1982) 



Selte 1 von 2 



internahonaler recherchenbericht 



Il Btlenales AMenzeielim 

PCT/EP 98/01595 



C.(Fort8etzung) ALS WESENTLICH ANGESEHENE UNTERLAGEN 



Katogorie' Bezetehnung der Ver6ffentlichur>g, aowert ertorderllch unter Angabe der In Betracht kommendGn Teile 



Betr. Anspruch Nr. 



VLADIMIROV V P, ALEKSANDROV E V: 
"Correction problem with a constraint on 
the number of Impulses" 
COSMIC RESEARCH, 

Bd. 18, Nr. 2, Marz 1980 - April 1980. 

Selten 128-137, XP002070356 

slehe Zusammenfassung 

slehe Selte 128, Absatz 1 - Absatz 2 

slehe Selte 134, Zelle 23 - Selte 135, 

Zelle 25 

DE 195 19 732 A (DAIMLER BENZ AEROSPACE 

AG) S.Dezember 1996 

In der Anmeldung erwShnt 

slehe das ganze Dokument 

US 5 335 179 A (BOKA JEFFREY B ET AL) 

2. August 1994 

slehe Zusammenfassung 

slehe Spalte 5, Zelle 30 - Zelle 66 

slehe Abblldungen 1,4 



1.2,10, 
11 



1.2 



US 5 140 525 A (SHANKAR UDAY J 

18. August 1992 

slehe Zusammenfassung 

BRONSTEIN I N, SEMENDADJEW K A: 
"Taschenbuch der Mathematik" 
1989 , VERLA6 HARRI DEUTSCH 
XP002070357 

slehe Selte 695 - Selte 717 



ET AL) 



1.2 



1-9,12 



THUN 



FomMW PCniSA/210 (FoilsaBung van eiaa 2) (Jlfl1992) 
BNSDOCID: ^V«> .9S4S058A1jL> 



Selte 2 von 2 



INTERNAnONALER RECHERCHENBERICHT 

Angaben zu Verdff entlio . .^n, die zur deiben Patent! amUie geh6ren 



Ir ttlonales Aktenzelchen 

PCT/EP 98/01595 



Im Recherchenbericht 
angefQhrtes Patentdokument 



Datum der 
Ver6ffentlichung 



Mitglted(er) der 
Patentfamilid 



Datum der 
Verdffentlichung 



us 


5130931 


A 


14-07-1992 


KEINE 




DE 


19519732 


A 


05-12-1996 


EP 0750239 A 


27-12-1996 


US 


5335179 


A 


02-08-1994 


KEINE 




US 


5140525 


A 


18-08-1992 


KEINE 





Fbtmblati PCTASAttl 0 <Anhan9 P&tanttamilteKJuB 19Se) 



aojtnneaA-f i «. 



